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Abstract 

In  this  paper,  the  integration  of  a  fuel  cell  into  the  power  system  is  treated  as  a  load  frequency  control  (LFC)  problem  with  the  fuel  cell 
acting  as  a  load  disturbance  source.  The  integration  of  a  fuel  cell  into  the  power  system  results  into  a  change  in  real  power.  But  changes  in 
real  power  affect  the  system  frequency.  Thus,  the  integration  will  result  into  a  change  of  frequency  of  the  synchronous  machines.  Hence,  we 
need  to  design  a  control  scheme  for  keeping  the  system  in  the  steady  state.  An  optimal  controller  based  on  the  disturbance  accommodation 
control  (DAC)  theory  is  proposed  for  this  load  frequency  control  problem.  For  demonstrating  the  effectiveness  of  the  proposed  controller, 
we  have  considered  a  two-area  power  system  with  the  fuel  cell  introduced  in  area  1 .  The  fuel  cell  is  considered  as  an  external  disturbance 
to  each  subsystem.  A  mathematical  model  is  derived  for  each  subsystem  and  based  upon  these  models  controllers  are  designed  for  keeping 
each  subsystem  stable,  which  in  turn  stabilizes  the  overall  system.  So,  the  proposed  controller  is  decentralized  in  nature.  To  account  for 
the  modeling  uncertainties,  an  observer  is  designed  to  estimate  each  subsystem’s  own  and  interfacing  variables.  The  controller  uses  these 
estimates  to  optimize  a  given  performance  index  and  allocate  generating  unit  outputs  according  to  the  requirements. 
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1.  Introduction 

The  increasing  demand  for  electric  power  forces  us  to 
consider  a  better  economic  alternative  rather  than  increasing 
the  load  on  transmission  lines.  The  option  of  constructing 
new  transmission  lines  is  restricted  because  in  spite  of  being 
feasible,  it  is  not  economically  sound. 

In  order  to  meet  the  local  demand,  an  economic  alternative 
in  the  form  of  small  generating  units  called  distributed  gen¬ 
erators  (DGs)  are  installed  on  the  power  system.  There  are 
different  types  of  distributed  generators  and  fuel  cell  (FC) 
is  one  of  them.  The  integration  of  a  fuel  cell  into  a  power 
system  will  cause  a  change  in  the  real  power  of  the  system. 
This  in  turn  affects  the  system  frequency  and  may  lead  to  in¬ 
stability  in  the  system.  Necessary  control  action  is  required 
to  bring  the  system  back  at  its  normal  operating  condition. 

In  this  paper,  a  new  type  of  optimal  controller  is  pre¬ 
sented  for  a  multi-area  power  system  using  some  results  from 
disturbance  accommodation  control  (DAC)  theory  [4,5].  A 


*  Corresponding  author.  Tel.:  +1-304-442-3205;  fax:  +1-304-442-3330. 
E-mail  address:  adavari@wvu.edu  (A.  Davari). 

0378-7753/$  -  see  front  matter  ©  2003  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2003.09.069 


power  system  controller  is  designed  to  integrate  a  polymer 
electrolyte  membrane  (PEM)  fuel  cell  into  the  power  sys¬ 
tem  using  the  disturbance-utilizing  mode  of  disturbance  ac¬ 
commodation  control.  This  paper  proposes  a  load  frequency 
control  (LFC)  scheme  using  disturbance  accommodation 
control. 

In  recent  years,  the  control  of  large  inter-connected  power 
systems  has  received  an  increasing  attention  [6,11].  It  has 
always  been  of  great  concern  to  the  power  industry  to  be  able 
to  maintain  electric  power  supply  in  the  face  of  unknown 
system  disturbance.  Load  frequency  control  or  automatic 
generation  control  (AGC)  is  the  mechanism  by  which  energy 
balance  is  satisfied  [3]. 

Under  normal  operating  conditions,  a  power  system  is 
continually  subjected  to  small  random  disturbances.  Typi¬ 
cal  examples  of  such  disturbances  are  small  changes  in  the 
scheduled  generation  of  one  machine  or  a  small  generator 
such  as  fuel  cell  added  to  the  network.  Integration  of  fuel  cell 
as  a  distributed  generator  to  the  power  system  can  be  con¬ 
sidered  as  a  disturbance.  The  transient  behavior  of  a  power 
system  following  a  disturbance  is  in  general  oscillatory.  If 
the  system  is  stable  then,  these  oscillations  damp  out  to¬ 
wards  the  normal  operating  point.  These  momentary  oscilla- 
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tions  yield  fluctuations  in  the  power  flow  over  transmission 
lines.  If  a  certain  tie  line  in  an  interconnected  power  system 
undergoes  excessive  power  fluctuations,  it  may  trip,  thereby 
disconnecting  the  two  groups  of  machines. 

The  main  objectives  of  the  research  presented  in  this  paper 
are: 

(i)  To  control  the  frequency  of  synchronous  machine  and 
tie  line  power  in  each  area  affected  due  to  real  power 
of  the  fuel  cell. 

(ii)  To  control  the  power  deviation  of  the  synchronous  ma¬ 
chine  in  each  area  affected  due  to  the  online  presence 
of  the  fuel  cell. 

This  paper  is  organized  as  follows.  Section  2  presents  a 
dynamic  model  of  the  power  system.  An  integrated  model 
of  the  power  system  with  a  fuel  cell  is  developed  in  Section 
3.  Section  4  provides  some  background  information  on  the 
disturbance  accommodation  control  theory.  In  Section  5  we 
design  an  optimal  controller  based  upon  the  disturbance  ac¬ 
commodation  control  theory  in  order  to  stabilize  the  inte¬ 
grated  power  system  with  the  fuel  cell  installed.  We  simulate 
the  integrated  power  system  with  the  controller  installed. 
The  results  obtained  from  the  simulation  are  presented  in 
Section  6.  Conclusions  are  discussed  in  Section  7. 


2.  Dynamic  model  of  the  power  system 

A  model  of  an  interconnected  power  system  consisting  of 
a  number  of  subsystems  or  control  areas  is  given  in  [2].  Each 
area  can  be  modeled  in  large  details  depending  on  the  gener¬ 
ator  models  and  their  prime  movers.  We  introduce  a  fuel  cell 
to  this  interconnected  power  system  as  depicted  in  Fig.  1. 

The  system  is  modeled  as  a  collection  of  independent  gen¬ 
eration  (GENCOs),  transmission  (TRANSCOs),  and  distri¬ 
bution  (DISCOs)  companies  where  the  fuel  cell  is  modeled 
as  a  power  generator.  GENCOs  produce  electric  power  that 
is  delivered  to  the  DISCOs  either  directly  or  through  the 
TRANSCOs.  In  the  structure  proposed,  the  DISCOs  are  to 
be  responsible  for  tracking  the  load  and  hence  performing 


load  frequency  control  task  by  securing  as  much  transmis¬ 
sion  and  generation  capacity  as  needed.  The  objective  of  the 
DISCO  is  to  supply  power  to  its  load  at  a  nominal  frequency. 
In  case  of  load  disturbance,  GENCO  1  will  adjust  its  out¬ 
put  accordingly  to  track  the  load  changes  and  maintain  the 
energy  balance. 

Using  the  linearized  dynamics  of  the  generators,  simplest 
models  of  speed  governors  and  turbines  associated  with  the 
generators,  the  power  flow  equations  from  generators  to 
distribution  system,  and  the  change  in  load  equations,  the 
two-area  test  system  are  modeled  as  follows. 

2.7.  Dynamics  of  the  generator 


Under  normal  operating  conditions,  the  relative  position 
of  the  rotor  axis  and  the  resultant  magnetic  field  axis  is  fixed 
[1,2].  The  angle  between  the  two  is  known  as  the  power 
angle  or  torque  angle.  During  any  disturbance,  the  rotor  will 
decelerate  or  accelerate  with  respect  to  the  synchronously 
rotating  air  gap,  and  a  relative  motion  begins. 

Applying  a  small  perturbation  to  the  swing  equation  of 
a  synchronous  machine,  the  dynamics  of  the  generator  are 
given  by 

2H\  d(A/j) 

- — -  =  A  JV,  —  A  Ei  —  D\  A/j 

jo  d  t 


2H2  d( A/2) 
fodt 


APm2  ~  AE2  -  P>2  A/2 


d(Afr) 
d  t 


2  7t  A  f 


where  A  represents  the  deviation  from  nominal  value,  77/  the 
constant  of  inertia,  79/  the  damping,/)  the  nominal  frequency, 
fi  the  actual  frequency,  8\  the  rotor  angle,  and  Pm  the  turbine 
(mechanical)  power. 


2.2.  Equations  for  speed  governors  and  turbines 


The  generators  are  equipped  with  speed  governors.  The 
simplest  model  of  a  speed  governor  and  turbine  associated 


Fig.  1.  Schematic  diagram  for  integration  of  a  fuel  cell  into  a  distribution  system. 
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with  the  generator  are  given  by 


d(A  Pyi) 
ch 

d(APMi) 
d  t 


A  Pyi  Khi  / 

Thl  Thi  V 

A  Pyii  +  Km  i 
Tui  Tut 


APrefi 

A  Pyi 


where  Py  is  the  steam  valve  power,  Tm  and  7h  the  time  con¬ 
stants  of  turbine  and  governor,  Km  and  the  gain  constants 
of  turbine  and  governor,  Rj  the  droop  characteristics,  and 
Prefi  the  reference  set-point. 


We  define 
Vi  =  \Vi\ZASi 
V2  =  |E2|ZA<$2 
Y\\  =  Gu  +jB\\ 
Yu  =  \Yu\^n 


The  above  equation  yields 

P\  =  |Vi|2Gii  +  |Vi||V2||yi2|cos(Aa1  -  A92-012) 

(6) 

Q  i  =  |Vi|2Bh  +  I  Vi  1 1 V21 1 3^12 1  sin(A3,  —  A32  —  012) 


2.3.  Power  angle  equations 


The  relation  between  rotor  dynamics  and  the  swing  equa¬ 
tion  is  expressed  by  [2]  as 


2H  d25 
&>s  d  t2 


In  the  swing  Eq.  (3)  for  the  generator,  the  input  mechanical 
power  from  the  prime  mover,  Pm ,  will  be  considered  con¬ 
stant.  This  is  a  reasonable  assumption  because  conditions 
in  the  electrical  network  can  be  expected  to  change  before 
the  control  governor  can  cause  the  turbine  to  react.  Since 
Pm  is  constant,  the  electrical  power  output  PQ  will  deter¬ 
mine  whether  the  rotor  accelerates,  decelerates,  or  remains 
at  synchronous  speed.  When  PQ  equals  Pm ,  the  machine  op¬ 
erates  at  steady- state  synchronous  speed;  when  Pe  changes 
from  this  value,  the  rotor  deviates  from  synchronous  speed. 
Conditions  on  the  transmission  and  distribution  networks 
and  the  loads  on  the  system  to  which  the  generator  supplies 
power  determine  changes  in  PQ.  Electrical  network  distur¬ 
bances  resulting  from  severe  load  changes,  network  faults, 
or  circuit-breaker  operations  may  cause  the  generator  output 
Pq  to  change  rapidly. 

Our  fundamental  assumption  in  this  work  is  that  the  effect 
of  machine  speed  variations  upon  the  generated  voltage  is 
negligible  so  that  the  manner  in  which  Pe  changes  are  deter¬ 
mined  by  the  power-flow  equations  applicable  to  the  state  of 
the  electrical  network  and  by  the  model  chosen  to  represent 
the  electrical  behavior  of  the  machine.  Each  synchronous 
machine  is  represented  for  transient  stability  studies  by  its 
transient  internal  voltage  Vi  in  series  with  the  transient  reac¬ 
tance  Xt,  Vt  is  the  terminal  voltage.  This  corresponds  to  the 
steady- state  representation  in  which  synchronous  reactance 
X  is  in  series  with  the  internal  or  no-load  voltage  Vi. 

Let  P  and  Q  represent  real  and  reactive  power  of  the 
system. 

We  have, 

N 

Pk  +jQk  =  VKJ2(YknVn)^<pin  +  A Sn  -  A S{  (4) 

n  —  1 


Let  us  assume  the  initial  condition  as  k  =  1  and  N  =  2. 
Therefore,  we  have 


We  assume 

A 8  =  ASi  -  A 82  (7) 


and  define  a  new  angle  y  such  that 

Tt 

Y  =  <t>  12  -  7 

The  above  equation  then  changes  to 

APe  =  APC  +  APmax  sin (5  —  y)  (8) 


The  above  equation  is  often  called  as  the  power  angle  equa¬ 
tion.  The  parameters  Pc,  PmSLX ,  and  y  are  constants  for  a  given 
network  configuration  and  constant  voltage  magnitudes  IVil, 
IV2I.  When  the  network  is  considered  without  resistance  all 
the  elements  of  Tbus  are  susceptances,  and  therefore  G\\  and 
y  are  zero.  Then,  the  power  angle  equation  becomes 


AP\  =  T\  sin(A<5i  — 


IVil  |V3|  cos(i5j 
Ai 


A53) 

-50) 


=  P 


max 


This  is  indeed  the  familiar  equation,  which  applies  for  pure 
reactance  network  given  by 

APe  =  AnaxSin(<$)  (10) 


where 


max 


Vj  \\Y2  \  cos  (3) 
Xi 


Similarly,  when  we  apply  the  power  angle  equation  to  bus 
2,  we  get 


A P2  =  T2  sin(A<$2  —  A<$3) 


(11) 


where 

IV2IIV3I  cos(SO-SO) 

1 2  =  - 

x2 

The  following  properties  hold  in  a  well  designed  and  prop¬ 
erly  operated  power  system: 

(i)  The  angular  differences  (5i  —  52)  between  typical  buses 
are  so  small  that 


Pi  +  JQ 1  =  V^YnVif  +  Vi(Ti2V2)* 


COS (5;  -  Sj)  =  1; 
sin (5,-  —  Sj)  =  (Si  -  Sj) 
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(ii)  The  line  susceptances  Bij  are  many  times  larger  than  2.4.  Model  for  the  power  system 
the  line  conductance  Gij  so  that 

p  .  „  /0  _  .  The  state  space  model  of  the  power  system  as  derived 

Gijsm(8i  -  8j)  ^  BijCos(8i  -  8j)  c.uu  •  •  u 

J  J  J  J  from  the  above  equations  is  given  by 

x  =  Ax  +  Bu  +  Gw 


xT  =  [A/i  APmi  A  Pyl  A8\  —  A82  A/2  APm2  A  Py2] 

U  =  A/^gfi 


wT  =  [A  PL] 


A  = 


where 
—  D  1 


a  = 


TP.  = 


Bt  = 


gt  = 


0 

—Khx 

RiThi 

2tx 

0 


0 

0 

T1T2 
T\  +  T2 

2  Hi 


fo 


TP\  Tpx 


0 

0 

0 

0 

0 


0 

K  M\ 

Tm\  Tm  1 

1 


1 


0  0 


Kff 

~Ti 


Th i 


0 

0 

0 

0 


—a 


Tp  ! 
0 


0 


0 


Tp2 

0 

0 


0 


0 


0 


—2tx 


-a  —  Z>2 


Tp , 


0 


~kh2 

BlT Hi 


L  (^i  +  r2)rPl 


000 


000 


Kh2 


Th 


2  j 


-r2 


0 

0 

0 

0 


1 


Tp , 


1 


0 


(14) 


0 

0 

0 

0 

0 

14  Mi 


Tm 9  Tm, 


1 


Tho 


(T\  +  r2)rPl 


0  0 


(iii)  The  reactive  power  Qi  injected  into  any  bus  I  of  the 
system  during  normal  operation  is  much  less  than  the 
reactive  power  which  would  flow  if  all  lines  from  that 
bus  were  short  circuited  to  reference,  i.e. 

Qi  l^l2% 


Therefore,  power  flows  between  GENCO  and  DISCO  are 
given  by 

AP^TiiASi- AS3) 

AP2  =  T2(AS2- AS3)  (12) 

A P2  =  -T2( ASi  -  A S2)  +  T2( ASi  -  A S3) 


where 

rji  _  |Vj||V3|  _ 

1  i  —  —  ^  max 

A/ 

The  change  in  load  is  expressed  by 
A PL  =  AP\  +  A P2  =  (T\  +  T2)(ASi  -  A S3) 

-T2(ASi  -  AS2) 

APl  T2 

(A Si  -  A53)  =  -  A52) 

7 1  T  f 2  7i  +  12 


(13) 


3.  Integrated  model  of  power  system  with  fuel  cell 


In  Section  2.4,  we  derived  a  state  variable  model  for  the 
power  system.  This  model  can  be  represented  in  the  gener¬ 
alized  form  as 


x  =  f(x,  u,  t) 
y  =  z(x,  u) 


(15) 


After  the  integration  of  the  fuel  cell  into  the  power  system, 
we  can  remodel  the  overall  system  as 


x  =  fix  +  Xp,  u ,  t) 
y  =  z(x-\-  x f,  u ) 


(16) 


where  xp  is  a  new  variable  introduced  as  a  result  of  the 
integration  of  the  fuel  cell. 

We  have  chosen  a  PEM  fuel  cell  for  its  easy  and  safe 
operational  modes,  low-temperature  gradient,  less  chance 
of  catalyst  poisoning  and  a  wide  scope  of  application  in 
power  distribution  system  [8-10].  A  detailed  modeling  and 
simulation  of  the  dynamic  behavior  of  this  fuel  cell  has  been 
carried  out  in  [10]. 
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The  result  of  the  chemical  reactions  inside  a  fuel  cell  is  the 
reversible  single  electrode  potential,  Er  given  by  the  Nernst 
equation  as 

?*)  In  (  P^° 

nF  J  \pu2  x  Jp^2 


where  E°  is  the  standard  electrode  potential,  R  the  gas  con¬ 
stant  (8.3144  J/molK),  T  the  temperature  in  Kelvin  scale,  n 
the  number  of  electrons  per  reacting  ion  or  molecule,  F  rep¬ 
resents  Faraday’s  constant  (96,500 C/mol),  pu2o  the  partial 
pressure  of  water,  pn2  the  partial  pressure  of  hydrogen,  and 
po2  the  partial  pressure  of  oxygen. 

A  chemical  polarization  of  cell  voltage  takes  place  due  to 
the  individual  chemical  reactions  occurring  at  the  electrodes 
of  the  fuel  cell.  This  is  given  by  the  Tafel  equation  as 


APchem 


RT 

anF 


In  i0  +  2.303 


RT 

anF 


log  (/) 


(18) 


where  a  is  the  activity  coefficient,  i0  the  exchange  current 
density,  and  I  the  applied  current  density. 

There  is  a  change  in  the  concentration  gradient  at  the  elec¬ 
trodes  due  to  the  continuous  chemical  reactions.  The  effects 
of  these  changes  are  expressed  in  the  form  of  concentration 
polarization  as  follows: 


AEC 


one 


RT 

nF 


In 


/l 


(19) 


The  total  power  output  from  the  fuel  cell  is  given  as 


P 1  cell  —  V0I0 


(24) 


where  V0  is  the  output  voltage  of  a  single-phase  inverter  and 
IQ  the  output  current. 

Using  the  voltage  dynamics  of  a  single-phase  bridge  in¬ 
verter  as  well  as  the  stack  current  dynamics  of  the  fuel  cell, 
we  get 


P 1  cell  = 


chem 


Vt-  A£t  /n,FAPT(t)  \ 
nRT  (dr  /  1 
,  ^  PT(t)  RTS  v\ 

+  [~l^-^imiXi-moXi) 

xVs(niF)cosa/IoRLin  /  /  >|  B 

nRTs  I  1  oinF  \  4) 


where 


(25) 


A  = 


B  = 


PORT  jn  /  pa 

pa  nF  V po 
/p  nFD 
II  -  /  A(1  -  a) 


In 


II 


where  7l  is  the  limiting  current  density  and  I  the  applied 
current  density. 

As  a  result  of  the  electrochemical  reactions,  there  is  a 
certain  amount  of  change  in  specific  conductivity,  leading 
to  an  additional  loss  in  the  potential  given  by 

nFD  //L-/\ 

A£Rest  =  ...  .  In  A—  (20) 

A(1  -  U)  \  II  J 

where  A  represents  the  equivalent  conductance  of  the  react¬ 
ing  ion  (cm2/ohm-equiv),  t[  the  transference  number,  and  D 
the  diffusion  constant. 

So,  the  sum  total  of  all  pertinent  components  of  polariza¬ 
tion  (APp)  occurring  in  a  single  cell  is  given  by 

A  A  —  AFchem  T  AE'conc  T  APrcsi  (21) 


Although  the  appearance  of  the  fuel  cell  is  only  restricted  to 
area  1,  there  may  be  a  possibility  that  the  effects  of  control 
actions  and/or  the  fuel  cell  itself  may  have  some  impact 
on  area  2  through  interconnections.  These  are  taken  into 
consideration  while  developing  a  model  for  the  integrated 
power  system.  We  develop  a  model  for  the  fuel  cell  with  the 
states  being  A/FC,  APichem,  APiceu,  APtie-FC,  zi,  zi,  Z3-  It 
is  very  important  to  note  that  zi,  Z2,  Z3  denote  the  impact  of 
the  fuel  cell  on  area  2. 

The  state  space  model  of  the  fuel  cell  is  given  as  follows 

[10]: 

Ff  =  Apvp  +  B^u 

(26) 

w  =  Cpvp 


Thus,  for  a  closed  circuit  cell  with  voltage  V,  operating  at  a 
load  current  /,  the  voltage  is  expressed  as  follows: 

V  =  (£r> A  -  A£t,a)  -  (Er,c  -  A£t,c)  -  i  J2  RJ 

V  =  £T,r  -  A£t,a  +  A£t,c  -  i  J2  Ri  (23) 

where  E^j  is  combined  standard  electrode  potential  of  anode 
and  cathode,  A£p,A  the  sum  total  of  all  pertinent  compo¬ 
nents  of  anodic  polarization  (a  positive  number),  and  AF/p,c 
the  sum  total  of  all  pertinent  components  of  cathode  polar¬ 
ization  (a  negative  number),  T  the  total  current  through  the 
cell,  and  Rj  represents  the  sum  of  all  internal  cell  resis¬ 
tance  including  electrolyte,  any  diaphragms,  and  resistances 
in  the  electrode  bodies. 


where 
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it.  This  state  space  model  is  given  as 


A  I 

B 

G 

Xi  = 

Xi  + 

u  T 

1 

& 

O 

_ 1 

By 

0 

3.1.  Linear  model  for  the  fuel  cell 


Xi  =  AjXi  +  Bju  +  Giw 


(28) 

(29) 


The  above  model  of  the  fuel  cell  is  nonlinear  and  cannot 
be  used  directly  for  controller  design  or  transient  stability 
analysis.  We  need  a  linear  model  for  the  fuel  cell  system. 
After  making  some  assumptions  based  on  practical  data  we 
transform  the  nonlinear  model  of  the  fuel  cell  into  a  linear 
one  given  as 


ip  =  Apvp  +  Byu 
w  =  Cpvp 


(27) 


where 
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3.2.  Integrated  system  model 

Using  relations  (14)  and  (27),  we  derive  a  state  space 
model  of  the  power  system  with  the  fuel  cell  integrated  in 


A\  and  B\  represent  the  integrated  power  system  model, 
A  and  B  the  original  power  system  model,  and  Ap  and 
By  the  fuel  cell  model  and  its  impacts  on  the  neighboring 
subsystem. 

The  model  represented  by  Eq.  (29)  will  be  used  for  de¬ 
signing  an  optimal  controller  for  a  smooth  integration  of  the 
fuel  cell. 


4.  Disturbance  accommodation  control  theory 

Design  of  optimal  controllers  based  upon  DAC  theory 
began  in  1967  [4,5]  and  has  since  been  developed  and  refined 
to  the  point  where  it  is  now  a  practical  and  general-purpose 
design  tool  ideally  suited  for  industrial  control  applications. 

Power  system  disturbances  are  not  totally  erratic,  but  have 
a  waveform  structure,  that  is  they  can  be  mathematically 
modeled  as  a  weighted  linear  combination  of  a  set  of  known 
basis  functions  given  as 

*f(0  —  Ci/i(0  +  Cifi(t)  +  C3/3W  +  •  •  •  (30) 

where  the  basis  functions  fit)  represent  the  various  wave¬ 
form  patterns. 

If  we  can  represent  the  unknown  input  disturbance  vp(0  by 
waveform  mode  characterization,  it  is  possible  to  disregard 
all  statistical  considerations,  random  process  theories  etc. 
and  proceed  to  design  a  physically  realizable  deterministic 
type  of  feedback  controller  based  upon  DAC  theory.  These 
types  of  controllers  are  remarkably  effective  in  coping  with 
the  specified  class  of  disturbance.  In  fact,  when  the  distur¬ 
bances  xp(0  have  waveform  structure,  DAC  typically  yields 
significantly  better  performance  than  the  so  called  stochas¬ 
tic  controllers  which  are  designed  by  considering  only  the 
experimentally  measured  long  term  statistical  properties  of 

*f(0  [5]. 

In  the  field  of  applied  control  technology,  it  is  a  well- 
known  fact  that  all  realistic  control  systems  operate  in  en¬ 
vironments  that  produce  system  disturbances  of  one  kind 
or  another.  Here,  the  term  disturbance  refers  to  those  spe- 
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cific  categories  of  the  system  inputs  that  are  not  accurately 
known  beforehand  and  which,  cannot  be  manipulated  by  the 
control  designer  i.e.  uncontrollable  inputs.  Disturbances  are 
an  important  factor  in  control  design  because  they  introduce 
unwanted  disruptions  to  the  otherwise  orderly  behavior  of 
the  controlled  system.  A  good  control  system  should  be  de¬ 
signed  in  such  a  way  that  it  maintains  the  given  control  spec¬ 
ifications  in  the  face  of  all  disturbances  that  might  act  on 
the  system  under  actual  operating  conditions. 

The  disturbance  accommodation  control  theory  allows  us 
to  systematically  design  multivariable  feedback  control  sys¬ 
tems  that  are  remarkably  effective  in  coping  with  both  tran¬ 
sient  and  persistent  disturbances  encountered  in  large  power 
systems. 

There  can  be  different  types  of  disturbance  accommoda¬ 
tion  controller  [4]. 

They  are: 

(1)  absorption  type, 

(2)  minimization  type,  and 

(3)  utilization  type. 

In  this  paper,  we  will  be  considering  the  utilization  type 
disturbance  accommodation  controller. 


5.  Control  of  power  distribution  system  using  the 
utilization  mode  of  DAC  theory 

The  main  idea  behind  the  working  of  a  disturbance  ac¬ 
commodation  controller  in  its  utilization  mode  is  that  distur¬ 
bances  might  sometimes  be  capable  of  producing  desirable 
effects  on  the  system  behavior.  From  a  conservative  point  of 
view  we  can  say  that  some  of  the  disturbances  can  be  con¬ 
structively  used  as  an  aid  in  carrying  out  the  primary  control 
objective.  For  example,  the  control  energy  and/or  the  tran¬ 
sition  time  required  to  bring  the  system  states  x(t)  to  some 
given  set  point  vsp  might  be  reduced  if  we  take  advantage 
of  the  natural  perturbed  motion  of  x(t)  caused  by  the  action 
of  disturbances.  Of  course,  we  have  to  manipulate  the  con¬ 
trol  u(t)  so  as  to  harness  and  exploit  the  useful  effects  in  the 
disturbances. 

If  the  disturbance  x^(t)  has  a  waveform  structure  that  can 
be  modeled,  we  can  employ  optimal  control  theory  to  design 


KX  =  (-A+BR 


Here,  maximum  utilization  of  the  fuel  cell  output  x^(t)  takes 
place.  We  regulate  the  system  frequency  deviation  A/j  and 
A/2  to  zero  while  keeping  the  time  interval  of  frequency 
deviation  as  small  as  possible. 

The  primary  control  objective  is  to  regulate  the  system 
output  y(t )  =  C(t)x(t)  from  any  initial  condition  y(t0)  = 
C(t0)x(t0)  to  the  given  set  point  y(T)  =  ysp  where  the  ter¬ 
mination  time  T  is  a  fixed  priori.  To  achieve  this  primary 
control  objective  and  simultaneously  make  maximum  uti¬ 
lization  of  the  fuel  cell  output  x?(t),  we  choose  the  control 
u(t)  to  minimize  the  quadratic  performance  criterion  given 
by  [4,5] 

1  1  [T 

J  =  -eft)  Se(t )  H —  /  [eft)  Qe{t)  +  uft )  Ru(t)]  d t 

2  2  A) 

(31) 


where  S ,  Q(t ),  and  R(t)  are  positive  definite  symmetric 
weighting  matrices  chosen  by  the  designer,  e(t)  denotes  the 
set  point  error  e(t)  =  ysp  —  y(t)  and  [t0  T]  is  the  specified 
time  interval. 

Now,  we  introduce  the  augmented  state  vector  x  = 
[x  ysp  vp  ]T-  So,  the  composite  set  with  the  condition 
ysp  =  0  can  be  written  as 


A  0  GCf 

”  B~ 
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The  set  point  error  e  can  be  expressed  in  terms  of  x  as 


e  = 


-CIO 


x 


(33) 


The  linear  quadratic  optimal  control  theory  is  applied  to  de¬ 
termine  the  control,  which  minimizes  (31).  The  final  expres¬ 
sion  for  the  disturbance  accommodating  controller  is 

u(x,xF,ysp,t)  =  -(R~lB'Kx)'x  -  (R~l  B'  Ky)ysp 

-(R~lB'Kx)x f  (34) 


where  Kx,  Ky,  Kz  are  time  varying  matrices  to  be  computed 
off  line  by  numerical  solution  of  the  set  of  matrix  differential 
equations  given  as 

1 B'  Kx)'  Kx  —  KXA  —  C'QC;  KX(T)  =  C'SC 


Ky  =  (—A  +  BR~l  B'  Kx)'  Ky  +  C'Q;  Ky(T)  =  -C'S  (35) 

K-  =  (—A  +  BR~l B' Kx)' Kz  -  KZAF  -  KXC¥ ;  KZ(T)  =  0 


a  disturbance  accommodating  controller  that  makes  optimal 
utilization  of  the  disturbances  [4]. 

5.1.  Integration  of  fuel  cell  into  power  distribution  system 
using  disturbance  accommodating  controller 

An  optimal  controller  is  designed  to  achieve  the  smooth 
integration  of  a  fuel  cell  into  the  power  distribution  system. 


A  part  of  the  control  signal  is  governed  by  the  Area  Control 
Error  which  is  generated  by  the  system  during  load  changes. 
This  is  represented  by  the  term  —  (R~lB'Kx)'x.  Also  part  of 
the  control  signal  is  governed  by  the  fuel  cell  output  since  we 
feedback  the  fuel  cell  output.  This  is  represented  by  the  term 
—  (R~lB'Kx)x p.  The  term  —  (R~l Bf Ky)ySp  governs  the  set 
point  associated  with  the  system. 

Solving  Eqs.  (34)  and  (35)  yields  [7], 
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u(x,  x¥,ysp,t)  = 


-495.94  -166.82  -10.59  -253.21  126.21  28.95  1.16 


y 


+ 


0  -101  2843  0  0  -10  0 


xp  243  ySp 


(36) 


The  variables  v  and  x¥  cannot  be  measured  directly  from  the 
system.  So  for  the  purpose  of  implementation  we  replace  v 
and  x¥  by  x  and  %,  which  are  the  outputs  of  an  observer. 

5.2.  Observer  design 


where  u  and  y  are,  respectively,  the  actual  output  and  the 
control  input  of  the  original  plant,  and  K\ ,  K 2  the  gain  ma¬ 
trices. 

These  gain  matrices  K\  and  K2  are  designed  by  examining 
the  error  dynamics  associated  with  the  observer.  We  define 


A  large-scale  power  system  consists  of  in  general,  a  num¬ 
ber  of  interconnected  subsystems.  Although  interconnection 
is  necessary  for  the  economical  operation  of  the  overall  sys¬ 
tem,  there  is  a  possibility  that  control  actions  and/or  exter¬ 
nal  disturbances  occurring  in  one  subsystem  will  be  trans¬ 
mitted  via  the  interconnections  to  the  neighboring  systems 
and  thereby  producing  undesirable  effects.  Hence,  in  an  in¬ 
terconnected  power  system,  not  all  the  state  variables  are 
measurable,  and  the  interface  variables  cannot  be  obtained 
from  local  measurements.  An  observer  is  used  to  overcome 
this  limitation  by  estimating  the  required  state  and  interface 
variables  using  only  available  measurements. 

The  output  of  the  observer  x  and  x¥  are  obtained  from 


€x  =  v  —  v 

*y  =  y  -  y 

Then,  the  error  dynamics  are  given  as 


A  +  KiC 

fc¥ 

_  f;v_ 

k2c 

D 

€y 

(38) 


(39) 


where  K\  and  K2  are  chosen  so  in  such  a  way  that  the 
transient  solution  of  the  above  equation  will  approach  zero. 

Applying  linear  observer  design  techniques  along  with 
pole  placement  methodology  to  the  above  equations,  the 


/V 

V 

5p 


A  +  KiC 
K2C 


(37) 


following  gain  matrices  are  obtained 
-57610;  6060;  -58000;  -155220;  -53710;  9830;  10580 

-160;  30;  -10;  -181600;  -28210;  -2390;  -100 


(40) 


Fig.  2.  Integrated  power  system  with  the  controller  installed. 
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A  schematic  diagram  of  the  integrated  system  with  the 
controller  installed  is  shown  in  Fig.  2.  It  depicts  a  two-area 
distribution  system  undergoing  load  frequency  control.  G1 
and  G2  are  the  synchronous  machines  (generators)  of  two 
different  areas.  The  fuel  cell  is  connected  in  area  1.  The 
difference  between  the  proposed  scheme  and  conventional 
load  frequency  control  is  that  instead  of  feeding  back  the 
tie  line  power  and  frequency  deviation  we  will  be  feeding 
back  the  tie  line  power,  frequency  deviation,  and  the  fuel 
cell  output  back  to  the  generators  [7].  The  results  of  this 
innovative  scheme  are  presented  in  Section  6. 

The  main  functions  of  the  proposed  controller  are  as  fol¬ 
lows: 

(1)  To  keep  the  system  frequency  approximately  at  the  nom¬ 
inal  value  (60  Hz). 

(2)  To  maintain  the  tie-line  power  flow  at  scheduled  value. 

(3)  To  regulate  the  generator  outputs  such  that  each  area 
absorbs  its  own  load  changes. 

In  the  next  section,  the  simulation  results  are  given  to 
show  how  well  the  controller  achieves  these  objectives. 

6.  Simulation  results 

6.1.  Tracking  performance  of  the  integrated  power  system 
for  different  types  of  load 

Figs.  3  and  4  represent  the  simulation  results  of  the  power 
system  incorporating  fuel  cell  and  a  tracking  load.  Fig.  3 


represents  the  tracking  of  the  system  to  a  sine  load  where 
as  Fig.  4  represents  the  tracking  of  the  system  to  a  square 
load.  The  output  of  the  system  is  tracking  the  input  very 
well  except  for  the  overshoot  at  the  start. 

6.2.  Power  deviation  and  frequency  deviation  using  a 
disturbance  accommodating  controller 

When  the  fuel  cell  is  connected  to  the  power  system,  there 
will  be  some  oscillations  in  the  output  of  the  generators 
which  should  be  dampened  out  quickly  to  avoid  any  dam¬ 
age  to  the  appliances  connected  to  the  system.  Ideally,  the 
power  flow  through  the  tie  line  is  zero.  But  during  a  distur¬ 
bance  power  starts  flowing  through  the  tie  line.  This  flow  of 
power  gradually  becomes  zero  as  the  system  becomes  sta¬ 
ble.  Similarly,  during  normal  conditions  the  power  system 
operates  at  a  particular  frequency.  But  at  the  time  of  a  dis¬ 
turbance  some  deviation  from  this  normal  frequency  takes 
place. 

Fig.  5  represents  the  power  deviation  and  the  tie  line  power 
flow  when  the  fuel  cell  is  incorporated  into  the  power  system 
and  there  is  no  controller  installed. 

From  the  above  response,  one  can  see  that  the  system 
oscillates  rapidly  when  the  fuel  cell  is  added  to  the  system 
and  takes  about  100  s  to  become  stable.  So,  we  incorporate 
a  disturbance  accommodating  controller  to  bring  the  system 
back  to  its  normal  condition  as  soon  as  possible. 

Fig.  6  represents  the  power  deviation  and  the  tie  line 
power  flow  when  the  disturbance  accommodating  controller 
is  present.  The  system  takes  about  7  s  to  become  stable. 


Fig.  3.  Tracking  behavior  of  the  integrated  power  system  for  a  sinusoidal  load. 
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Fig.  4.  Tracking  behavior  of  the  integrated  power  system  for  a  square  load. 
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Fig.  5.  Power  deviation  and  tie  line  power  flow  without  the  disturbance  accommodating  controller. 


228 


A.  Paradkar  et  al.  /  Journal  of  Power  Sources  128  (2004)  218-230 


Fig.  6.  Power  deviation  and  tie  line  power  flow  with  the  disturbance  accommodating  controller  installed. 


So  we  see  that  the  system  returns  to  its  normal  operating 
condition  much  faster  than  the  previous  case  when  the  con¬ 
troller  was  not  installed.  This  implies  that  the  disturbance 
accommodating  controller  is  very  effective  in  controlling  the 
system. 


Fig.  7  represents  the  frequency  deviation  when  the  fuel 
cell  is  incorporated  into  the  power  distribution  system  with¬ 
out  any  controller.  As  real  power  from  the  fuel  cell  is  added 
to  the  system,  the  system  frequency  gradually  increases  and 
the  system  becomes  unstable. 


Fig.  7.  Frequency  deviation  without  the  disturbance  accommodating  controller. 


A.  Paradkar  et  al.  /  Journal  of  Power  Sources  128  (2004)  218-230 


229 


x  10'3 


Time 

Fig.  8.  Frequency  deviation  with  the  disturbance  accommodating  controller  installed. 


Fig.  8  represents  the  frequency  deviation  when  the  fuel 
cell  is  integrated  into  the  power  distribution  system  in  the 
presence  of  a  disturbance  accommodating  controller.  We  see 
that  the  frequency  deviation  for  both  machines  becomes  zero 
in  a  very  short  period  of  time. 

This  controller  achieves  an  efficient  response  with  sub¬ 
stantially  smaller  frequency  deviation.  The  designed  con¬ 
troller  is  robust  to  external  and  neighboring  disturbances  like 
small  changes  in  the  scheduled  generation  of  one  machine 
or  a  small  load  added  to  the  system  as  well  [7] . 


7.  Conclusions 

The  integration  of  a  fuel  cell  into  the  power  distribution 
system  is  achieved  in  a  systematic  way  by  applying  dis¬ 
turbance  accommodation  control  theory.  The  designed  con¬ 
troller  makes  optimal  constructive  use  of  the  fuel  cell  output 
to  enhance  the  overall  performance  of  the  integrated  sys¬ 
tem.  The  performance  of  the  controller  is  assessed  through 
simulation  of  a  two-area  power  system.  The  main  objective 
of  this  control  scheme  is  to  keep  the  power  system  stable 
when  a  fuel  cell  is  connected.  The  controller  achieves  a  fast 
response  with  substantially  small  frequency  deviation  while 
maintaining  the  energy  balance  between  generation  and  con¬ 
sumption. 


We  observe  that  the  overall  system  stability  depends  not 
only  on  the  dynamics  of  the  power  system  but  also  on  the 
dynamics  of  the  fuel  cell.  The  dynamics  of  the  fuel  cell 
does  not  affect  the  system  stability  when  it  is  offline,  but  it 
definitely  affects  the  overall  system  stability  when  it  is  on. 

The  controller  developed  in  this  paper  will  be  used  to 
carry  out  a  detailed  analysis  of  the  different  penetration  lev¬ 
els  of  the  fuel  cell  and  its  effects  on  the  power  quality  of 
the  overall  system.  As  the  penetration  level  of  the  fuel  cell 
increases,  the  power  quality  will  deteriorate  and  necessary 
control  action  must  be  taken.  According  to  existing  laws, 
the  supplied  power  must  meet  certain  standard  requirements 
regarding  harmonics  and  power  factor. 

The  developed  controller  will  also  be  very  useful  for  re¬ 
liability  analysis  of  the  overall  system.  We  also  need  to  de¬ 
termine  the  fault  tolerance  of  the  integrated  power  system 
during  large  disturbances. 

All  these  efforts  will  eventually  lead  to  the  development 
of  a  cost-effective,  decentralized  power  system  involving 
significant  penetration  of  distributed  energy  resources  like 
fuel  cells. 
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